ON AN ISOMONODROMY DEFORMATION 
EQUATION WITHOUT THE PAINLEVE PROPERTY 



BORIS DUBROVIN AND ANDREI KAPAEV 



Abstract. We show that the fourth order nonlinear ODE which 
controls the pole dynamics in the general solution of equation Pj 
compatible with the KdV equation exhibits two remarkable prop- 
, erties: 1) it governs the isomonodromy deformations of a 2 x 2 

matrix linear ODE with polynomial coefficients, and 2) it does 
not possesses the Painleve property. We also study the properties 
' of the Riemann-Hilbert problem associated to this ODE and find 

, its large t asymptotic solution for the physically interesting initial 

r <P ■ data. 
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1. Introduction 

The study of relationships between the theory of isomonodromic de- 
formations and the theory of differential equations satisfying the so- 
called Painleve property is a well-established branch of the analytic 
theory of differential equations in complex domain (see below a brief 
summary of the most important results obtained in this direction). One 
j> \ of the outputs of the present paper suggests that the above mentioned 

relationship is less straightforward than it was traditionally believed. 
We illustrate the point with an example of a fourth order ODE for a 
r-». ! function a = a(t) 

O ' a"" + 120(a')V - 120a' a" t - fV) 2 - f aa" + f> t = 0. (Pj 2 ' 1} ) 

m ■ 

This equation appeared in [22] in the study of pole loci of solutions to a 

(2 1) 

degenerate Gamier system. Our observation is that eq. PJ ' governs 
isomonodromic deformation of certain linear differential operator with 
polynomial coefficients (see eq. (II. 5p below). However, this equation 
does not satisfy the Painleve property as its general solution has third 
order branch points of the form 

a(t) = a - (t - 6) 1/3 + O ({t - 6) 5/3 ) , kC (1.1) 

(cf. [12]) where the location of the branch points depends on the choice 
of the solution (the so-called movable critical singularities) . 

In spite of such a somewhat surprising phenomenon the method of 
isomonodromic deformations proves to be almost as powerful in the 

I (2 1) I 

study of solutions to the equation Pj as in the case of classical 
Painleve equations. Namely, it is possible to derive the large t asymp- 
totics of solutions and, moreover, to describe the branching locus of a 

l 
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given solution in terms of a kind of a spectral problem for a quintic 
anharmonic oscillator. 



The equation Pj 2,1 ^ is of interest on its own. Namely, it describes the 



behavior of poles of solutions to another fourth order ODE 

Uxxxx + 10w 2 . + 20uu xx + 40 (u 3 - Qtu + Qx) = (1.2) 

usually denoted as P 2 . It is the second member of the so-called Pi hier- 
archy. The coefficients of this equation depend on t as on a parameter. 
It is well known (see, e.g., [30J) that the equation (ll.2p is compatible 
with the Korteweg - de Vries (KdV) dynamics 

u t + uu x + j^u xxx = 0. (1.3) 

The Laurent expansion of solutions u = u(x,t) to the system (II. 2p . 
(11.31) near a pole x = a(t) has the form 

1 



u 



(x-a(t)) 



2 



O (Or - a{t)f) (1.4) 



where the function a(t) solved eq. Pj 2 ' 1 ^ In this way one associates 
a multivalued solution a(t) to eq. Pj 2,1 ^ with any solution to eq. P 2 . 



The branch points (II. ip of a(t) correspond to the triple collisions of the 
poles of solutions to the P 2 . Thus, the above mentioned isomonodromy 



(2 1) 

realization of eq. P\ ' provides one with a tool for studying the KdV 
dynamics of poles of solutions to the P 2 equation. 

Particular solutions to the P 2 equation are of interest. We will con- 
centrate our attention on one of them, namely, the one that has no 
poles on the real axis x e IR (see [3], [TT], [12], [5] about importance of 
this special solution to the eq. P 2 ). Such a solution u(x,t) exists [7] 
for any real t and it is uniquely determined by its asymptotic behavior 
for large \x\. Using the developed techniques along with the isomon- 
odromic description [26J of the special solution to P 2 , we arrive at the 
asymptotic description for large t of poles of this special solution. 

Before proceeding to the formulation of main results let us first 
briefly recall the basic notions of the present paper, namely the Painleve 
property and the isomonodromic deformations. In the description of 
the historical framework we will manly follow the paper 



1.1. Painleve property. In 1866, L. Fuchs p.4] has shown that all the 
singular points of solutions to a linear ODE are among the singularities 
of its coefficients and thus are independent from the initial conditions. 
In the nonlinear case, the reasonable problem is to look for ODEs 
defining the families of functions, called the general solutions, which 
can be meromorphically extended to the universal covering space of 
a punctured Riemann surface with the punctures determined by the 
equation. In other words, the problem is to find ODEs whose general 
solutions are free from the branch points and essential singularities 



1 The connection of eq. Pj' ' with the KdV equation was not considered in 
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depending on the specific choice of the initial data. This property is 
called now the Painleve property, or the analytic Painleve property, and 
is obviously shared not only by the linear ODEs but also by the ODEs 
for elliptic functions. 

In [T5], L. Fuchs started the classification of the first order poly- 
nomial in u and v! ODEs, F(x, u, v!) = 0, with single-valued in x 
coefficients with respect to the Painleve property. H. Poincare [31] and 
P. Painleve [32J accomplished the analysis not finding new transcen- 
dental functions. 

In [33], P. Painleve revisited the L. Fuchs' idea extending the program 
of classification to the second order ODEs of the form u" = F(x, u, u') 
with F meromorphic in x and rational in u and u'. In the course of clas- 
sification of the 2nd order 1st degree ODEs modulo the Mobius trans- 
formation [TBI [21], it was found 50 equations which pass the Painleve 
a-test and thus called now the Painleve-Gambier equations. It occurs 
that all these equations can be either integrated in terms of the clas- 
sical linear transcendents, or elliptic functions, or reduce to one of the 
six exceptional classical Painleve equations Pi— Pyi- 

1.2. Isomonodromic deformations. The monodromy group to a 
linear ODE was first considered by B. Riemann [37], H. Schwarz [ID] 
and H. Poincare [35]. Apparently, it was L. Fuchs in [16J who for the 
first time set the problem of the deformations of the coefficients in a 
linear equation that leave the monodromy group unchanged. Namely, 
assuming that solutions of a linear ODE depend on an additional vari- 
able, he obtained a system of the first order PDEs the solutions must 
satisfy. 

More modern treatment of the isomonodromy deformations was de- 
veloped by R. Fuchs in |T7] • He has shown that the monodromy group 
of a scalar linear ODE with four Fuchsian singularities at A = 0, 1, oo, x 
and an apparent singular point at A = u does not depend on the lo- 
cation x of the fourth Fuchsian singular point if the location of the 
apparent singularity u depends on x according to a nonlinear 2nd or- 
der ODE. Later, the Fuchs' isomonodromic deformation equation was 
included as the sixth Painleve equation Pyi into the list of the classical 
Painleve equations. 

In 1912, L. Schlesinger [39] generalized the Fuchs' approach finding 
the equations of the isomonodromic deformations for arbitrary linear 
Fuchsian ODEs, while R. Gamier [19] presented the scalar second order 
linear ODEs with irregular singular points whose isomonodromic defor- 
mations are governed by the lower classical Painleve equations Pi-Py. 
In 1980, M. Jimbo, T. Miwa and K. Ueno [25J extended the theory of 
Gamier to the linear ODEs with generic irregular singularities. 

1.3. Painleve property of the isomonodromy deformation equa- 
tions. It is interesting that to the date of the mentioned above achieve- 
ments, the fact that these equations do indeed possess the Painleve 
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property was not proved rigorously even for the case of the classical 
Painleve equations. Thus the term "Painleve equation" loosely refers 
to various equations among which we mention the higher order ODEs in 
hierarchies associated with the classical Painleve equations, the higher 
order ODEs in the classifications by F. Bureau [2] and C. Cosgrove 
[8] based on the use of some Painleve tests, as well as the differen- 
tial, difference, g-difference and elliptic-difference equations found in 
the course of the study of the symmetries and geometry of the classical 
Painleve equations, see e.g. [31] and [38] . 

A general elegant approach to the Painleve property of the equations 
of the isomonodromic deformations was presented by T. Miwa [29] and 
B. Malgrange [28] in the early 80s of the last century (recall that the 
direct proof of the Painleve property of the classical Painleve equations 
appeared even later, see [22]). The approach by T. Miwa and B. Mal- 
grange is based on the use of the zero curvature representation and the 
Riemann-Hilbert correspondence. In fact, they have proved the an- 
alytic Painleve property of the isomonodromic deformation equations 
for arbitrary linear Fuchsian ODEs and for equations with unbranched 
irregular singular points. In this respect we also mention the papers 
of M. Inaba and M. Saito [23] who developed an algebro-geometric ap- 
proach to the geometric Painleve property of the isomonodromic defor- 
mations of the logarithmic and unramified irregular connections which 
also implies their analytic Painleve property. 



1.4. Plan of the paper and the main results. One can expect 
that the isomonodromy deformation equations for arbitrary rational 
connections possess the Painleve property. However, the naive induc- 
tion does not work. Namely, and this is our first result, w e show that 
the polynomial in all variables 4th order 1st degree ODE (PJ ' ) gov- 
erns the isomonodromic deformations of the following linear differential 
equation with polynomial coefficients 



dZ 



20 U 



^A 4 + a 3 A 3 + a 2 X 2 + ctiA + a 



Z 



30 



20 " 



;i.5) 



Here a = a(t), otk = oif.(t), k — 0, 1, 2, 3 are some smooth functions. 
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Proposition 1.1. The monodromy^ of the system (11.51) does not de- 
pend on t if and only if the coefficients have the form 



I 3 

a 3 = — —a! + «3, a 2 = y^ a ' — 603a' — t + a\ 

9 /3 , ry7„.0-./2 



ct x = + 27a°a' - 3[a° + 60(a°) 2 - t]a + a - 40a°t + a° 

a = - JL"' - ^ a /4 + 432«°a' 3 + 18 [t - 330(a°) 2 - a°] a' 2 
+3 [a - 130ajjt + 10800(a§) 3 + 90ajja!j + a?] a - 30a§ a 



+2W0(a° 3 ) 2 t - 54000(a 3 1 ) 4 - 900ajj(aj}) 2 - 30a?a° 



where al, a 2 , a® are some constants, together with the following ODE 
for the function a = a(t) 



a IV 



120a' 3 - 3600a°a' 2 + 120 (a° + 270(a°) 2 - t) a' 



40 

- — (a - lOOajjt + 6300(a£) 3 + 60a£a£ + a?) 



„ 200 , 2 200 

a a H c 

3 9 



+^ a ° a ' _ 6000(a°) 2 - — a° 2 = 0. (1.6) 
3 9 

After a change 

t^t + a° 2 - 30(ajj) 2 , a ^ a + 10a° 3 t - a? + 40a° 2 a° 3 - 300(a£) 3 



(2 1) 

the equation (I1.6P reduces to the P\ ' equation. 



(21)1 

So, our claim is that the equation P\ [ describes the isomonodromic 
deformations of the system (11.51) but, however, it does not possess the 
Painleve property as being satisfied by a 4-parameter Puiseux series in 
powers of [t — 6) 1 / 3 . Nevertheless, as we show, this equation can be 
effectively analyzed using the isomonodromy deformation techniques 
developed during last decades. In particular, the ramification points 
t = b (see eq. (II. ip above) can be determined by a kind of a spectral 
problem for the quintic anharmonic oscillator 

(1.7) 

U(X) = ^A 5 - oA 3 + a A 2 + (*f b 2 + § a 9 )A - f* a b + fo u . 

Here ao, ag, an are the coefficients of the Puiseaux expansion of a(t) 
(see eq. (14.11) below). Namely, the coefficients of the quintic polynomial 
must be chosen in such a way that the equation (I1.7P possesses solutions 
exponentially decaying on certain contours in the complex plane, look 
for details in Remark 14.11 



2 As the system (|1.5[) has only an irregular singularity at infinity, the isomon- 
odromicity means independence from t of the Stokes multipliers of the system - see 
below. 
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In Section 2 we discuss the second member of the Pi hierarchy, equa- 
tion Pj, as an equation that describes the isomonodromic solutions to 
the KdV equation, and find a nonlinear ODE, equation Pj , that 
controls its pole dynamics. In Section 3, we derive the above isomon- 

(2 1) 

odromic representation for the eq. P\ ' presenting a regularization 
of the linear system for Pj along its singularity locus. In Section 4, 

(2 1) 

we consider the Puiseux series solution for P\ ' and the singularity 
reduction of the corresponding linear system at the branch points. 
In section 5, we set the Riemann-Hilbert problems for all the previ- 
ously introduced wave functions and discuss existence and uniqueness 
of their solutions. In Section 6, we present the asymptotic analysis of 
the Riemann-Hilbert problems corresponding to a physically interest- 
ing special solution of Pj and p[ 2,1 ' ) implementing the steepest-descent 
method introduced by Deift and Zhou. Our main result in this section 
is the description of the large t asymptotics of the singularity locus as 
a theta-divisor on a modulated elliptic curve, 

3x + §£i(zo) = lt~ 7/i {{n + \)uJ a {xo) + (m + \)u b (x )), 

where n G Z, m G Z + and xq = a(t)t~ 3 ^ 2 , the function £i(a?o) is de- 
termined by equations ( I6.6p - (16.8I) . u a (xo) and oj^xq) are the period 
integrals ( 16. lip on the elliptic curve (jfi.iup . In Section 7, we find ex- 
plicitly the pole dynamics of the special solution in a vicinity of the 
attracting point Xq = |^|, 

+ «(«+i)^L + Or 7 / 2 ln 2 t), (1.8) 

where m G Z + and n G Z enumerate the points of the pole lattice. 
Thus, asymptotically for large t — >■ +oo, the poles of the special solution 
to Pj never collide. It would be interesting to study their collisions for 
finite time. In Section 8, we briefly discuss our results and various open 
problems. 

2. Linear system system for Pj 2 and its singularity locus 

2.1. Linear system for Pj. The 4th order 1st degree polynomial 
ODE, 

u xxxx + I0u 2 x + 20uu xx + 40(u 3 - 6tu + 6x) = 0, (2.1) 

usually denoted as Pj , is the second member of the so-called Pi hierar- 
chy. This equation governs the isomonodromy deformations of a linear 
polynomial ODE described by the system 

®x = AV, ^ = 5*, ^ t = C7#, (2.2) 
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where the connection matrices are explicitly written using the genera- 
tors of jeI(2,C), a 3 = a+ = (§J) and a_ = (°°), 

A = 2io( [ _4m ^ ~ {^uu x + u xxx )]a 3 

+ \8X 2 + 8uX + (12u 2 + 2u xx - 120t)lff+ 
+ [8A 3 - 8uX 2 - (Au 2 + 2u xx + 120*) A 

+ (16m 3 - 2u 2 x + 4uu xx + 240x)] a_ } , 
B =a + + (X-2u)a-, (2.3b) 
C =lu x a 3 - |(A + - |(A 2 -uX- 2u 2 - Im^)^-, (2.3c) 
Besides ( 12. ip which is the compatibility condition 

[d x -A,d x -B] = 
of (I2.3al) and (I2.3bp . the system (12. 3p also implies the KdV equation, 

u t + uu x + ±u xxx = 0, (2.4) 

and equations which follow from (12.11) and (12. 4ft . The solutions of Pj 
(12.11) compatible with the KdV equation (12.4D are called the isomon- 
odromic Pj solutions to KdV. 

2.2. Laurent series solutions to Pj compatible with the KdV 
equation. It is known that equation Pj passes Painleve tests as being 
presented in the list by C. Cosgrove [S] under the symbol F-V with 
y = —u and the parameters a = 240£, k = 240 and (3 = 0. See [UJ 
for a proof of the Painleve property based on the Riemann-Hilbert 
correspondence. 

Below, we are especially interested in the 4-parameter series solution 
to Pj with the following initial terms (because Pj is polynomial in 
all its variables, the construction of the complete formal series via a 
recurrence relation is straightforward), 



"(■'•) : --, ^7 + ^ c k (x - a) k , 

[x — a) z / — ' 



c 1 = 0, c 2 = 3(c 2 -2t), c 4 = f a - 10c 3 + ±f%t, 

C 5 = 3 — |c C 3 , C 7 = y(tC 3 — C ), 

ci, Co, C3, Cg are arbitrary. (2.5) 

This 4-parameter series is compatible with the KdV equation ( 12. 4p 
if the coefficients cq, C3, cq depend on t in a particular way described by 
the function a(t), 

c = a, c 3 = 2a", c 6 = -\a"' - 3(a') 4 + 12(a') 2 t - 12t 2 , ( )' = 4- 

6 at 

(2.6) 

In its turn, the pole position a(t) must satisfy the equation Pj , 

+ 120(a') 3 a" - U0a'a"t - ^f(a') 2 - faa" + 2g°t = 0. (2.7) 
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Besides the pole dynamics, the function x = a(t) parametrizes the 
singularity locus for the linear system (12. 2p . (12. 3p . 

Remark 2.1. Equation P 2 admits also the 3-parameter series solution, 

3 °° 

u{x) = — — 9 + y~] c k {x - a) k , 

(x — a) z ' 

v ' k=2 

C 2 = — |t, C 3 = 0, C 4 = — |j, C 5 = —1, C 7 = 0, Cg = — jgjt, 

? — A(-^0_n 2 _L_ 11? + _i_ J6_+3\ ~ _ 20 ~ 
° 10 13 I 1323" " r 7 C6t " r 343 6 / ' Cn ~~ 441 

a,c 6 ,c 8 are arbitrary, (2.8) 

which, however, is not compatible with the KdV equation (12. 4p . Clearly 
it corresponds to the triple collisions of the poles (12. 5p . 



3. Singularity reduction in the linear system for Pj and 
the isomonodromy property of pj 2 '^ 

3.1. Singularity reduction. It is remarkable, that the linear system 
for P 2 can be regularized along the singularity locus x = a(t): 

Theorem 3.1. Let u(x) be the J^-parameter Laurent series solution 
Ii2.5\) for equation P 2 . Then the gauge transformation Z = R^fr with 
the gauge matrix 

R(\) = 2^ 3 (A - 2c )^(x - a )-^ 3 j={a 3 + ajx 
x (A - 2c )-^2^(x - a)-° 3 ±(a 3 + a x ){x - a)-^\ a x = 

(3.1) 

regularizes the connections d\ — A and d x — B along the singularity 
locus x = a(t) and, if 

Co = a' and C3 = 2a" , 

it regularizes also the connection dt — C along the same locus. 

Proof. The proof is straightforward by inspection. □ 

Remark 3.1. The above gauge transformation is constructed as a se- 
quence of the so-called shearing transformations. Namely, substituting 
the initial terms of the relevant Laurent series instead of u into the 
matrix A, one finds the leading order singular term proportional to 
the sum —(x — a)~ 4 a_ + (2 — a)~ 3 cr 3 + (x — a)~ 2 a + . Conjugation by 
(x — a) _CT3//2 "shares" the singularity order among the off-diagonal en- 
tries. The matrix coefficient of the term of order (D((x — a)~ 3 ) becomes 
nilpotent, and its conjugation with ^75(03 + o~\) turns it to the normal 
cr_-form. This brings into play the lower order singularities, and the 
procedure can be repeated. 
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Remark 3.2. It is known that the singularity reductions exist also for 
the linear systems associated with the classical Painleve equations. 
However the construction of the corresponding gauge transformation 
does not rely on a sequence of the shearing transformations. In con- 
trast, the required gauge transformations simply turn the 1st order 
rank 2 differential systems into scalar 2nd order linear ODEs. 

Conjecture 3.1. The singularity reduction exists for an arbitrary isomon- 
odromy system. 

(2 1) 

3.2. Isomonodromy representation for PJ' . The linear system 
for Pj regularized along its singularity locus x = a(t) yields the linear 
system for equation Pj 2 ' 1 ^. 

(2 1) 

Theorem 3.2. Equation P\ ' controls the isomonodromy deforma- 
tions of a linear matrix ODE with the polynomial coefficients described 
by the system 

Z x = AZ, Z t = CZ, (3.2) 
with the coefficient matrices 

^ = ~li a//cr 3 + (^j^ 4 + a 3^ 3 + a^A 2 + «iA + a )cr + + (^A - a 3 )cr_, 
a 3 = -j^a\ a 2 = ^(a') 2 ~t, oti = -^j(a') 3 + 3a't + a, 
ao = 3aa ' _ 27( a /)4 + |t( a ') 2 + fee + ft 2 , 
C = [-§A 3 + ftA 2 + ftA + ft]a + - 

ft = 2a', ft = lOt - 9(a') 2 , ft = 36(a') 3 - GOa't - 10a. (3.3) 

Proof. By inspection, the compatibility condition, A t — C\ + [A, C] = 0, 
yields 

c 6 = -la'" - 3(a') 4 + 12(a')H - 12t 2 , 

and 

a"" + 120(a') V - 120a'a"t - ^(a') 2 - faa" + ^ft = 0, 
i.e. equations for c 6 in (I2.6P and Pj ' (12.71) . □ 



4. PUISEUX SERIES FOR a{t) AND THE SINGULARITY REDUCTION 
OF THE LINEAR SYSTEM AT THE BRANCH POINT OF a(t) 

(2 1) 

It is easy to verify that equation Pj ' admits the formal series solu- 
tion in powers of (t — 6) 1 / 3 with the following initial terms (again, since 
the equation is polynomial in all variables, the recurrence relation for 
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its coefficients is straightforward), 

oo 
k=0 

a\ = —1, a2 = a 3 = 0,4 = 0, 05 = — 1&, a6 = 0, 

a 7 = ^ao) &8 — — |) &10 = 0, 012 = — -Jg-&> • • • , 

b,a ,ag,au are arbitrary. (4.1) 

The system (13 .2L H3. 31) is singular, however regularizable, at the 
branch point t = b of a(t): 

Theorem 4.1. Let a(t) be the J^-parameter Puiseux series solution 
d^.l| ) for equation Pp. Then the rational in A and t gauge trans- 
formation X = QZ with the gauge matrix 

Q(A) = a 3 (±fb - X 2 Y^\t - b)^^ 2 (a 3 + a,)x 

x 2- CT3 A-^3(i|o 6 _ x 2 y a3 (t - b)^ 3 ^(a 3 + (Ti)x 

regularizes the connection d\ — A at the branch point t = b. 
Here 01, cr 2 , a 3 are Pauli matrices, 

ai= (l J)' a2= (z "0)' as= (o -!)■ 

The new function X(X) satisfies the linear matrix ODE equivalent 
to the anharmonic oscillator equation (II. 7p . which does not admit any 
continuous isomonodromic deformation, 

Xx X" 1 — I ? 30 1 
AaA - \U(X) J ' 

/7(A) = ^A 5 - bX 3 + a A 2 + ( 3 -§ b 2 + § a 9 )A - a 6 + fo u . (4.3) 

Remark 4.1. The anharmonic oscillator equation (11. 7ft can be regarded 
as a stationary Schrodinger equation. Thus it is possible to introduce 
a scattering problem and to map the set of the parameters b, ao, ag, 
an to the set of the scattering data. Indeed, let us introduce the Jost 

(k) 

solutions if)± uniquely determined by the asymptotics 

^ - A-*"£ ( ± \) e-, = ^ - |tA 3 / 2 + xA 1 ^, 

A ^00, argA = tt + ^ A;, fc = 0, ±1, ±2, ±3. 

Here, we define A 1 / 2 on the A complex plane cut along the positive part 
of the real line and choose its principal branch. Then the incident wave 
tp ~ ip^ produces several reflected and transmitted waves, 

ijj = t k ^ ] + r k ^\ argA = 7r+ 2 ffc, k = 0, ±1, . . . , ±3, 
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where the coefficients if. and are called the transmission and reflec- 
tion coefficients, respectively. Only 4 of the 14 coefficients are indepen- 
dent. Indeed, there exist six relations of the form r^ + i = r^, = t^-i, 
k = —2,0,2, which come from the preservation of the amplitude of 
the dominant solution in the sectors arg A G [tt + ^-k, n + ^-(k + 1)], 
k = —3, • • • , 2; the equality to = 1 normalizes the amplitude of the in- 
cident wave; two conditions £3 = r_3 = mean the solution is subdom- 
inant as A — > +00; the last condition r 3 = it_ 3 means the continuity of 
this subdominant solution across the positive part of the real line, 

h = r-3 = 0, t = *-i = 1, t 2 = h(= 0), 
n = r (= i), r_i = r_ 2 (= 0), r 3 = r 2 = it_ 2 = it- 3 {= i). (4.4) 

The four independent scattering coefficients, say, t±, r±, t_ 2 , r -2, can be 
used to reconstruct the parameters b, a , a 9 and a u determining the 
potential U(X). 

In (14.41) . the numbers in parentheses specify the values of the free 
transmission and reflection coefficients corresponding to the triple col- 
lisions of the poles of the special solution to Pj of our interest, see 
below. Along the oscillatory directions, the wave function correspond- 
ing to this special potential is given by the Jost solutions ifj^ (k = 

-1, -2, -3), iip { + ] (k = 1, 2, 3) and by the sum ^i 0) + iipf as A < 0. 
This implies another simple characterization of the potential of our in- 
terest. Defining A 1 / 2 on the plane cut along the negative part of the real 
axis, the Schrodinger equation ip" = ^U(X)if^ has the solution with the 
uniform in the sector arg A G (— 7r,7r) asymptotics ip — A~ 5 / 4 e -61 . Ex- 
istence of a solution with such asymptotics is tantamount to vanishing 
of four Stokes multipliers s±± = s± 2 = (see below). 

Below, however, we characterize X(\) adopting the Riemann-Hilbert 
problem instead of the scattering problem. 

Remark 4.2. Observe the successive reduction of the number of the 
deformation parameters the A-equations depend on: for the connection 
d\ — A, the space of the deformation parameters (t, x) is 2-dimensional, 
for d\ — A, the deformation parameter space reduces to the 1-dimensi- 
onal space of t, and for d\ — 3^0+ — Z7<7_, it is 0-dimensional. 

Remark 4.3. Besides the continuous isomonodromic deformations con- 
sidered above, it is possible to introduce discrete Poincare-like isomon- 
odromy mappings from one to another branch of a(t) and from one to 
another point of the lattice (b,a(b)). 

Remark 4.4. It is natural to denote the above listed connections by 
the symbols l| 2,2 \ Lj 1 and Lj 2 ' '', respectively, to reflect explicitly 
the number of the continuous deformation parameters involved. The 
equations of the isomonodromic deformations are thus denoted by the 
symbols p| 2,2 ^ (same as P 2 ), P^'^and Pj 2,0 ^(the latter equation is a 
Poincare-like mapping of the lattice (b,a(b))). 



12 BORIS DUBROVIN AND ANDREI KAPAEV 

5. RlEMANN— HlLBERT PROBLEMS FOR L,f' 2) , L^'^AND L| 2 ' 0) 

The boundary Riemann-Hilbert (RH) problem consists in finding a 
piece-wise holomorphic function by its prescribed analytic properties: 

i) asymptotics at some marked points; 

ii) discontinuity properties across a piece-wisely oriented graph. 

5.1. Asymptotics of the canonical solutions ty k (X), Z k {X) and 
X k (X). Each of the above linear matrix ODEs for Z and X has 
one irregular singularity at A = oo and no other singular points. We 
introduce the formal series solutions that represent the asymptotic be- 
havior of the genuine solutions to these linear ODEs in the interior of 
particular sectors near A = oo. 

Since all the canonical asymptotics differ from each other in a ratio- 
nal left diagonal multiplier, it is convenient to unify the notations and 
distinguish all three cases introducing a parameter v: 

*M(A) = A^(a 3 + a x ){I + 0{\-^))e^\ 
fl= JLA^-ftA^ + xA 1 / 2 , 

A^oo, argA G (-f + ^fc, f + fk), k G Z, 
ve {-1,3,-5}, $(- 1 ) = ^, $( 3 ) = Z, $ ( ~ 5) =X (5.1) 

Here, the principal branch of the square root of A is chosen. 

Remark 5.1. Using ( 15. ip for the asymptotic solutions Z k and X k , one 
has to take into account the reductions x = a(t) for v — 3 and t = b, x = 
a(b) for v = —5. Below however, using ( 15. ip to set the corresponding 
RH problems, we assume that the exponential 9 depends on arbitrary 
complex deformation parameters t and x because it is not allowed to 
use the unknown values a(t) and b, a(b) as the data in the RH problems 
formulated below. As the result, it is necessary to remember that the 
solvability domains for the RH problems shrink to the lines x = a(t) 
for v = 3 and to a lattice (t, x) = (b, a(b)) for v = —5. 

5.2. Stokes multipliers. The canonical solutions differ from each other 
in a triangular right matrix multipliers called the Stokes matrices, 

$ fc+1 (A) = MX)S k , Sh-i = (I > S ^ = Q 2k J) > ( 5 - 2 ) 

moreover 

Sk+7 = s k , s k + s k+2 + s k s k+1 s k+2 = -i{l + s k+4 s k+5 ), k G Z. (5.3) 

Remark 5.2. Since the gauge transformations R(X) and Q(X) are ra- 
tional, they do not affect the Stokes matrices. 

Remark 5.3. The Stokes multipliers s k are the first integrals of equa- 
tions Pj 2,2) (J22) and Pj 2,1) (12T7D . 
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The solutions of the RH problems below are constructed using col- 
lections of the canonical solutions in the above presented sectors in a 
way to have the uniform asymptotics in a vicinity of infinity. The nor- 
malization of the RH problem however is rather subtle thing. To this 
aim, observe the prolongated canonical asymptotics 

±{a 3 + a 1 )X^ k (X)e- e ^ = (I + § a, + 0(X^ 2 ))e x - 1/2d ^ , (5.4) 

= (' - l>i + M>') 2 - + <?(a- 5/2 ))x 

x expiX-^d^s + X- 3/2 d 3 a 3 + X- 2 dJ}, (5.5) 

and 

±(a 3 + a 1 )X^X k (X)e- 9 ^ = (/ + - ^ + 0{X- 7 ' 2 ))* 

x exp{X- 1/2 d l( r 3 + X- 3/2 d 3 a 3 + X^' 2 d b a 3 }. (5.6) 

In the above right hand side exponential factors, the parameters dj can 
be expressed in terms of the coefficients of the relevant equation. For 
instance, d\ in (15 .4p is one of two Hamiltonians of Pj. However, these 
particular relations are not important for us at this stage. 
Now, we have the following 

Riemann— Hilbert problem 1. Given the complex values of the pa- 
rameters x, t and Sk, k e Z, satisfying Ii5.3\) . and denoting 

0=^X 7 l 2 -\tX 3 / 2 + xX 1 ' 2 , (5.7) 

find the piece-wise holomorphic 2x2 matrix function $^(A) ; v E 
{ — 1, 3, —5}, with the following properties: 

(1) 

lim X(a 3 + a 1 )X-^ s <l> {u \X)e- e(73 = I, (5.8) 

moreover 

• there exist a constant d\ such that 

^(a 3 + a l )X-^^\X)e- 9 ^ = 

= I + X- 1 ' 2 d 1 a 3 + 0(X- l ) ] (5.9) 

• if v 6 {3, —5}, then there exist constants di, d 3 such that 

-±(a 3 + a 1 )A-5 ff3 2^ 3 $M(A)e- (e+dlA_1/2)<T3 = 

= I + X- l Cl + X- 3 / 2 d 3 a 3 + 0(X- 2 ) (5.10) 
with some constant matrix c\; 
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• if v = —5, then there exist constants di,d 3 and d 5 such 
that 

^(a 3 + a 1 )A-t^2^3 $ H (A)e -(^ 1 A-v 2+rf 3A-3/ 2 ) (J 3 = 

= / + \- 2 c 2 + \- 5/2 d 5 a 3 + (D(\- 3 ) (5.11) 

with some constant matrix c 2 ; 

(2) ||$ H ( A )II < const as A ^ 0; 

(3) on the union of eight rays 7 = p U (u 7 k=1 ~fk-4) , where jk — 
{X E C: argA = ^fk}, k = -3, -2, . . . , 2, 3, and p = {A e 
C: argA = 7r} ; all oriented towards infinity, the jump condition 
holds true, 

$^(A) = $ H (A)S(A), (5.12) 



where $V ( A ) and ®- ( A ) are limits of <& {u) {\) on 7 from the 
left and from the right, respectively, and the piece-wise constant 
matrix 5(A) is given by equations 




FIGURE 5.1. The jump contour 7 for the RH problem [T] 
and canonical solutions $j(A), j = —3, —2, . . . , 3, 4. 

5.2.1. Uniqueness of the solution to the RH problem \J\ The solution 
of the RH problem [1] is unique. Indeed, the scalar function det $^(A) 
is continuous across the set of rays 7 and is bounded at the origin. 
Therefore, det$ M (A) is an entire function, and taking into account 
(15. 8 p and applying the Liouville theorem, det $(^(A) = —1. Assume for 
the moment that there are two solutions of the RH problem HJ $^(A) 
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and $M(A). Consider their ratio, x(A) = = $M(A)($ M (A)) -1 . 
As it is easy to see, x(A) is continuous across all the rays of the set 7 
and is bounded at the origin, thus x(A) is an entire function. Using 

63, 

ltaA-^ X (A)Ai«=lim( c ; /2 W ; /2 ) =/ ' 

Since all the entries a, 6, c, d are entire functions in A, the Liouville 
theorem yields the ambiguity 

$ (l/) (A) i-> 4> (l/) (A) = P M (A)$ (l/) (A), 

p(-D(A) = / + c a_, P< 3 >(A) = / + ( Co A + Cl )a+, 

p(- 5 )(A) = / + ( Co A 2 + ciA + c 2 K, (5.14) 

where Cj are arbitrary constants. These ambiguities however are elim- 
inated using the asymptotic conditions in (15.91) . (15.1 Op and (15. lip . 

5.3. Solvability of the RH problems [I] and the Malgrange di- 
visor. Given jump matrices, the set of points (x,t) at which the RH 
problem is not solvable is called the Malgrange divisor. It coincides 
with the zero locus of a holomorphic r-function by Miwa and with the 
singularity locus of the isomonodromy deformation equation. 

For what follows, it is convenient to articulate our assumptions on 
the singularities and critical points of equations Pj and P[ 2,1 ' ) . 

(2 2) 

Conjecture 5.1. 1) Equation Pj ' has no other movable singularities 
except for the movable poles (12. 5p satisfying (12.61) . (12. 7p . and their triple 
collisions (12. 8p . 2) Equation Pj ' has no other movable singularities 
or critical points except for the branch points (14. ip . 

These assumptions mean that the smooth branches of the Malgrange 
divisor for the RH problem [TJ with v — — 1 (corresponding to Pj ) are 
parameterized by equations x = a(t) and thus coincide with the solv- 

(2 1) 

ability set for the RH problem [T] with v = 3 (corresponding to Pj ). 
Similarly, the vertices of the Malgrange divisor of the RH problem with 
v = — 1 correspond to the branch points t = b of a(t) and therefore 
coincide with the solvability set of the RH problem with v = —5. 

In other words, Con] ect ure 15.11 implies the following remarkable prop- 
erties of the domains of solvability of all three RH problems: 

(1) these domains do not pairwise intersect; 

(2) all these domains together cover the deformation parameter 
space C x C 3 (t, x). 

5.4. Normalized non-homogeneous and homogeneous RH prob- 
lems. In this section, we describe another interesting feature of the 
family of three RH problems. Roughly speaking, any solution of the 
non-homogeneous RH problem with the bigger value of \u\ allows one to 
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x=a 3 (t) 



Figure 5.2. The scheme of the Malgrange divisor and 
the domains of the RH problems solvability for v = 
-1,3,-5. 

construct infinitely many solutions of the homogeneous problems with 
the smaller value of \u\. 

Precise formulation of this property requires some more accuracy. 
First of all, let us introduce the RH problems equivalent to those above 
but unbranched and normalized at infinity. To this aim, write 

A = C 4 , ^(0 = ^3 + aOC^^CV^, (5-15) 
and define the piece-wise holomorphic functions, 

& V) {Q = *P(Q, argCe(£(*-l),5fc), ^ = -1,3,-5. (5.16) 

These functions solve the non-homogeneous RH problems on the union 
of rays Ul=_ 13 £k, £k = {( E C: argC = jjk} with a singular point at 
the origin and normalized to the unit at infinity. The corresponding ho- 
mogeneous RH problem differs from the non-homogeneous counterpart 
in the asymptotics of that vanishes as £ — > oo. 

Slightly abusing our notations, we formulate the non-homogeneous 
and homogeneous RH problems as follows: 

Riemann— Hilbert problem 2. Find a piece-wise holomorphic func- 
tion $P(£) ( resp., $o (C))j v £ { — 1)3,— 5} ; with the following prop- 
erties: 

(1) lim^oo $^(C) = I (resp., lim^ <lj%) = 0); 

(2) across the rays £k = {( E C : argC = jjk}, k = — 13, . . . , 13, 14, 
all oriented towards infinity, the following jump conditions hold, 

&:\() = ^\0e e ^S k e- e ^, Ce4, f? = ^C 14 -K 6 + < 2 ; 

(5.17) 

(3) ||C^^((73 + tri)$M(C)||<cofiat, as ( -> 0. 
The following theorem holds true 
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Theorem 5.1. Any solution (Q to the non-homogeneous RH prob- 
lem [H with v = 3 yields infinitely many solutions $g ^ (£) to £/ie ho- 
mogeneous RH problem^ with v — — 1. ylny solution &\ (C) to the 
non-homogeneous RH problem [H z/ = —5 yields infinitely many 
solutions $o^(C) ^° ^ e homogeneous RH problem^ with v = 3, — 1. 

Proof. We prove the first part of the theorem. The proof of the second 
part is similar. 

Consider a solution to the RH problem [2] for v = 3. It has 

the asymptotics 

$5 3) (O = / + 0(r 2 ), c^oo, 

^ 3) (C) = ^(^3 + ^i)r 3,T3 (Mo + 0(C 4 )), C^o, (eu , 

where M Q is a constant matrix, det M 7^ 0, and Wo = {( G C: arg£ G 
(0,^)}. Then 



cti), p,q — const, 



has the same jump properties as $j (£) an d the asymptotics of the 
form, 

^- 1> (0 = o(r 2 ), c^oo, 

*t\o = ^(^i)r(o g ^(M +o(c 4 )), c^o, (6 Wo . 

Thus $0 (C) is the solution of the homogeneous RH problem [2] with 
u=-l. □ 

6. Large t asymptotics of a special solution of equations 

P 2 AND Pf^ 

In this section, we construct large t asymptotic solution to the RH 
problem [1] for v — — 1 and ^ = 3 corresponding to a special solution 
of equation P 2 . This special solution u(x,t) has the asymptotics it ~ 
=F-y/ 6|x| as x —>■ ±00 and is real and regular on the real line for any 
tel. The physical importance of this solution for t — from the point 
of view of string theory was justified in [3J. In [TlT [T2"] . this solution 
appeared in the study of the problem of universality of critical behavior 
of solutions to Hamiltonian perturbations of hyperbolic PDEs. 

The above mentioned properties uniquely distinguish this special so- 
lution. In [26], the characterization of this solution in terms of the 
Stokes multipliers of the associated linear system was found, 



S-2 = S-l = Si = s 2 = 0, s_ 3 = s = s 3 = -i. (6.1) 
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6.1. Large t asymptotic spectral curve for the special solution 

to Pj. Assume that t > is large and x G C is such that 

Xo ■= x t~ 3/2 = 0(1) as t ->• +oo. (6.2) 

Our starting point is the large t asymptotics of the spectral curve, 
det( / u - r 5/4 A(A)) = 0, where A = t 1/2 £, i.e. 

5 

S = M 5 ~ s^ 3 + M^ 2 + j- D^ + j- D =<k - 60- (6.3) 



k=l 

)2 



The asymptotic analysis of various degenerated solutions of Pj per- 
formed in (20], shows that the topological properties of the asymptotic 
spectral curve are significantly different in the interior of the domains 
D± and D , see Figure IBTTl For x G D±, the spectral curve has genus 0, 

^ = ^(e-&) 2 (e-^) 2 (£-60, (6.4) 

where the double branch points and the simple branch point 
satisfy the conditions, 



6,, = -\£ k ± V15 - fel e k - 24£ fc + 48x = 0, 

and the ambiguity in choice of the root of the cubic equation for the 
simple branch point is fixed demanding that, for xq G K, 

sgn(xo) Re ^ J (j,(z) dz^j > 0, xq — > ±oo. 

The conditions above are consistent with the quasi-stationary asymp- 
totic behavior of the special solution to P 2 ( 12. ip . 

u(x, t) ~ t 1//2 vo, ^0 — 6i>o + Qx = 0, 

Vq ~ —^6x0, x — > ±oo, (6.5) 
where the real on the real line branch of \/6xn is chosen. 




Figure 6.1. The discriminant set for the special solu- 
tion of equation P 2 . In the interior of the domains D + 
and D_, the asymptotics of the special solution has genus 
g = 0, and in the domain D , it has genus g = 1. 
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Let us indicate some interesting points in the Xq complex plane. 

At the point xq = — 2\/3 = D$ n D- n R, two double branch points 
of the asymptotic spectral curve coalesce, so that £1,2,3,4 = — V% and 
£ 5 = 4\/3- Generically, at the asymptotically quadruple branch point, 
the local solution of the RH problem can be approximated using the 
Gamier- Jimbo-Miwa \l/ function for the second Painleve transcendent 
Pii; [EE [25] . For the RH problem we are studying, the relevant local 
approximate solution corresponds to a Hastings-McLeod solution to 
Pii ®. 

At any x G D + including its boundary, the spectral curve has one 
simple and two double branch points. The precise asymptotic location 
of these branch points can be found for xq = = D H D + PI R, where 

£ 3 = -^p £1,2 = and £ 4j5 = At this point, as well as at 

any point of _D + including its boundary, the leading order asymptotic 
solution of the RH problem can be expressed in elementary functions, 
cf. [1] for the real line case. 

Observe also the point Xq = 2\/3 G D + where two double branch 
points coalesce. The asymptotic branch points corresponding to this 
point are £3 = — 4a/3 and £1,2,4,5 = v3- Generically, the quadruple 
degeneration corresponds to the appearance of Pn. However, for the 
RH problem we consider here, the relevant Painleve function is trivial, 
and the asymptotic solution to the RH problem remains elementary. 

In the interior part of the domain D Q , the large t asymptotic solution 
to the RH problem is constructed on the model elliptic curve, 

m 2 = 9m(£ - ^i) 2 (£ - m - m - &), (e.e) 

where the branch points j = 3, 4, 5, are determined by the values of 
£1 and x as the roots of the cubic equation 

e + 2^ 2 + (3£? - 30)£ + 4^ - 60£x + 30a;o = 0. (6.7) 

The double branch point £1 is determined as a function of x G -Do by 
the system of Boutroux equations, see , 

Re / fjt(z) dz = 0, Re / fi(z) dz = 0, (6.8) 

supplemented by the described above boundary conditions on D±C\D . 

Remark 6.1. As x approaches the real segment (— 2\/3, |^), all as- 
ymptotic branch points become real and satisfy the inequalities £3 < 
£1 < £4 < £5. Thus the second of the equations (16.81) trivializes while 
the first of these equations turns into the condition fi(z) dz = 
obtained in [SB] in the analysis of the Whitham equations and used in 
[1] to study the same special solution on the real line. 

6.2. Steepest-descent analysis of the RH problem. The strategy 
of the steepest-descent asymptotic analysis by Deift and Zhou P, \TU\ 
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of the RH problem involves several standard steps: 1) transforma- 
tion of the jump graph to the steepest-descent directions of a suitable 
g-function; 2) construction of the local approximate solutions (para- 
metrices); 3) matching all the parametrices into a global parametrix; 
4) proof that the global parametrix indeed approximates the genuine 
solution to the original RH problem. 

Since all the above mentioned steps are well explained in the litera- 
ture (see e.g. [13]), below, we omit unnecessary details. 

6.2.1. Transformation of the jump graph to the steepest-descent direc- 
tions. According the steepest-descent strategy, we first transform the 
jump contour for each of the RH problems to the steepest-descent graph 
for the exponential exp{ J fi(z) dz}, see Figure [6721 




Figure 6.2. Typical anti-Stokes lines for the special so- 
lution of Pj as Xq e -D , Imx > 0. 

Observe that, in the special case (16.11) . s± 2 = s±i = 0, s = s± 3 = —i, 
the jump graph depicted on Figure 15.11 can be transformed to the one 
shown in Figure 16.31 where 

1, J) , 2+ := (J 7) * <-> 

6.2.2. Model elliptic curve and abelian integrals. The large t asymp- 
totics of \I/ ( A) corresponding to the Stokes multipliers (16.11) is con- 
structed on the Riemann surface T of the model elliptic curve 

^ 2 = (£-6)(£-£4)(£-60, (6.io) 

glued of two copies of the complex £-plane cut along [£ 5 , £4] U [£3, —00), 
see Figure 16741 
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Figure 6.3. The jump graph for the RH problem in the 
special case s± 2 = s±i = 0, s = s± 3 = —i. 




Figure 6.4. The Riemann surface V and the basis of 
cycles a, b. 



Define the complete elliptic integrals 

A,B=<f M£) £iM£) df, 

«/ a, ft •/ a, ft 



W a .ft 



and abelian integrals 
9(0 = 



r = — , imr > (J, 



/i(z)dz, £/(£) 



j_ /" € dz 

"a 7 ?5 tw(z) 



(6.11) 



(6.12) 



We define the integral g(£) on the upper sheet of the Riemann surface 
T cut along the sum of intervals [£5, £4] U [£4, £3] U [£3, —00). 
The Boutroux conditions (16. 8ft imply 



A,B E iR. 

Observe the following properties of <?(£) and £/(£): 



(6.13) 
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(1) as £ — > oo, 

U(0 = U oo + O(C 1/2 ), C/oo = -|r, 

(6.14) 

(2) g(£) and t/(£) are discontinuous across the broken line [£5, £4] U 
[^4^3] U [£3, —00) oriented from £ 5 to infinity, moreover 

£g(£ 5 ,£ 4 ): g+(0+9-(0 = 0, U + {0 + U40 = 0, 
£ G (£4,6): <7+(£) -</-(£) = -A, C7+(0 - MO = -i, 
£ G (6, 00) : <?+(£) + #_(£) = —B, M£) + MO = -r. 

(6.15) 

6.3. "External" parametrix. In this subsection, following [13], we 
solve the permutation RH problem on the segments [£5, £4] U [£3, —00) 
whose solution gives a leading order contribution to the solution of the 
above RH problem. 

Riemann Hilbert problem 3. Given t ^> 1 and the complex value 
of the parameter x , find a piece-wise holomorphic 2x2 matrix function 
<&^(£) with the following properties: 

(1) lim -L(<t 3 + a 1 )C^-t CT 3$M(£) e -* 7/4 ^3 = j (6 16) 

^ = ^ 7/2 -K 3/2 + ^£ 1/2 ; (6.17) 
across the union of segments (£5, £4) U (£3, —00) oriented as indicated, 
the jump condition holds true, 

e&M(fl = *W(fl«ri, £ G (£ 5 ,£ 4 )u(£ 3 ,-oo) (6.18) 

where $^(£) and $ H (£) are limits o/$^(£) on £/ie segments from 
the left and from the right, respectively, see Figure W7c\ 

We do not impose any conditions on the behavior of $ ( - l/ ^(£) at the 
points £ = £-,-, j = 3,4,5. As the result, the solution of the model 
RH problem [3] is determined up to a left rational matrix multiplier 
with possible poles at £ = £j, j = 3,4,5, and certain asymptotics at 
infinity. Below, we will use this rational multiplier to prove or disprove 
the asymptotic solvability of the original RH problem [TJ 

6.3.1. (3 -factor. Consider the principal branches of the functions A,(£), 

/MO = (£ - 6r 1/4 (£ - £ 4 ) 1/4 (£ - £ 5 r 1/4 (6.19) 
&(£) = (£ - £ 3 ) 1/4 (£ - £ 4 ) 1/4 (£ - £ 5 ) 1/4 (6.20) 

both defined on the £ complex plane with cuts along the broken line 
[£5, £4]U[£4£3]U[£3, —00). These functions solve the scalar RH problems: 

(1) &(£) = £i(l + 0(£- 1 )), as £^00; 
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(2) the discontinuity of /3j,(£) across the oriented contour [£5, £4] U 
[^3] U [£ 3 , —00) is described by the conditions 

£e(£ 5 ,£ 4 ): P± 1 (£) = i0I 1 (Z), #(0 = -iPsiO, 

£e(£ 4 ,£ 3 ): ^(0 = ^1(0, /#(£) = -& _ (£), 

£e (£ 3 ,-oc): (3^(0 = 1^), /#(£) = ^(O- (6.21) 

6.3.2. The Riemann theta function and the Baker- Akhiezer functions. 
Define one more function, 

K(0 = t 7/ \g(0-9oo) + s u (u(0-u oo ), goo = -\B, U 0O = -\t, 

(6.22) 

where the parameter S u is defined by 

£_i = -t 7 ^A, 5 3 = -t 7 ' A A + in. (6.23) 

This function /i„(£) has the following obvious properties, 

(1) as £ — > +00, 

h v (Z) = t 7 "0 + O(C 1,2 y, (6.24) 

(2) h v (£) is discontinuous across the broken line [£5, £4] U [£4, £3] U 
[£ 3 , —00) oriented from £ 5 to infinity, moreover 

£e(£ 5 ,£ 4 ): /*+(£) + / i ;(£) = * 7/4 £ + ^, 
£e(£ 4 ,£ 3 ): K(Z)-K(Z) = -t 7/A A-6 v , 

£G(£ 3 ,oo): fc+(0 + MO = 0- (6-25) 

Using the Riemann theta-function, 0(z) = ^ n e™ 2 ^ 2 ™^ define 
the matrix function $l Sj4 ^ (£ ) , 



X 



/ e(t/(g)+y+<M //.x e(-t/(Q+v+^) \ 

\ e(t/(0) ©(-^(0) c 2lsjy 

Here the parameters V, 4> u and the factors c,-(£), c*(£) are defined by 

^ = -^ ?/4 M-S); 
0-i = ^, 03 = 0; 



c 2 (£) = /?I 1 2 (£) i ^y, ^(£) = /T 2 (£)|^fy if ^ + ^ = 5, 

n,mGZ. (6.27) 
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It can be shown that det & BA) (£) = -1, the function flBT26D 
satisfies (I6.16P and (I6.18P and thus is one of the solutions of the RH 
problem [31 Any other solution to this RH problem has the form of the 
product 

¥J 3A \0 = R^H BA) (0, (6-28) 
where R u {0 is rational with poles at £ = £j, j = 3,4,5, and satisfies 
the asymptotic condition R u (£,) = / + (9(^ _1 ) as £ — > oo. 

6.4. Local RH problem solution near the branch points £ = 

£j, j = 5,4,3. As it is well known, near a single branch points, the 
relevant boundary- value problem can be solved using the classical Airy 
functions, see e.g. [T3] . 

6.4.1. RH problem for the Airy functions. Define the Wronsky matrix 
of the Airy functions [I], 

where 

Vl {z) = Ai (z), v 2 (z) = e i2n/3 M (e l2lT/3 z). (6.30) 
Besides Z (z), introduce auxiliary functions 

Z. x (z) = Z (z)(S + )- 1 , Z x {z) = Z Q (z)Z_, Z 2 (z) = Z 1 {z)E + , 

E + = (J T) , - J) ^ (6.3D 

By construction pQ, 

Z y («) = ^ CT3/4 ^(^3 + <n)(/ + 0(z-t))eI 23/2CT3 , 
2^00, z£w 3 -{zeC: arg2 G (-tt + f j, § + fj)}- (6.32) 
Assemble the piece-wise holomorphic functions Z^'(z), j = 5,4,3, 



Z ( %) 



'Z_ 1 (^)G i , argze(-7r,-f) 
Zq^G,, argzG (-f ,0), 



Z 1 (z)Gj, argz G (0, 
K Z 2 (z)Gj, argze 



'2tt 

- 3 ■ 



3 

3 / 
7rl 



(6.33) 



where 

G5 = G 3 = /, G 4 = za 3 . 
Observing that the jump matrices of Z^'(z) coincide with those 
across the lines emanating from the node points j = 5,4,3, in 
Figure 16.31 we are ready to construct the relevant local parametrices, 

¥j\0 = B^(0Z {j) (z {j) (0), J =5,4,3, ^ = -1,3. (6.34) 

Here Bj v) (0, 

j = 5, 4, 3, are holomorphic in some finite neighborhoods 
of £ = £j matrices, and 2 = are changes of variables biholomor- 

phic in some finite neighborhoods of £ = 
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6.4.2. Determination of the local change z = zV'(£). This biholomor- 
phic change of variables has to be chosen in a way to ensure that the 
global parametrix, see below, has small enough jumps as t — > +00, 
namely it must satisfy the condition 

|(* ( %)) 3/2 = t 7 ^(g(0-9m+Su(U(0-Um+oW, le-Cil = const. 

The biholomorphicity condition is satisfied with the choice 

/ft \ 2/3 

= t V«(* ^ ft z)dz j i j = 5j4j3) (6.35) 

where /f(f) has the form (16. 3p . (16. 6 p with the branch points fj, 

The unique abelian differential holomorphic on the elliptic curve w 2 = 

(£ -&)(£- is 



Thus the 1-parameter deformation of the model degenerated curve that 
respects the degeneration is generated by the vector field 

& = (6-36) 



dD u> a \ dbi Si db ) ' 

Finally, we find the elliptic curve satisfying the asymptotic conditions, 

mo + t- 7/A ^ = m + o{t-w{t - orn 

as |f - £i| > dr 7/8 , |f - fj| > (7,-r 7 / 4 , Cj = const, j = 3, 4, 5, 

oj a , b = u a , b + C(r 7 / 2 ), f = r + o(r 7 / 2 ), 

^ + r 7 /% = A + G(r 7 / 2 ) , s + t^ 7/4 r^ = B + G(r 7 / 2 ) . (6.37) 

6.5. Global parametrix. The global approximate solution to the RH 
problem for 3>^(f ), v = —1, 3, is a piece-wise analytic matrix function 
^i/(0 defined as follows, 



*,(0 



*? A) (0 = ^(0$ (BA) (0, Cj = EC: If - 0| < r}, 

0<r<imin|6-e i |, fc,j = 3,4,5, (6.38) 

see Figure 1631 

The exact solution is constructed using the correction function x(f), 

x(0 = *(£)* _1 (0, (6.39) 
that satisfies the following RH problem: 
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Figure 6.5. The construction of the global parametrix. 

(1) the limit 

exists and is diagonal; 

(2) across the contour 7 shown in Figure 16.61 the jump condition 
holds, = X-(0 H (0 where 

/f(0 = *-(0S(0*+(0- 




FIGURE 6.6. The jump contour 7 for the correction 
function %(£)■ 

To apply the L 2 -theory to the latter RH problem [H], all jump ma- 
trices H(£) across the jump contour 7 in Figure 16.61 must satisfy the 
estimate — ||L 2 ( 7 )nc(7) — o(l) as t — > 00. For the infinite tails 

emanating from the circles, this fact holds true because the relevant 
jumps are uniformly exponentially small. 

The jumps across the circles centered at the branch points £j, j = 
3, 4, 5, can be made small as t is large if one adjusts the rational matrix 
Rv(0> v = — 1;3, and the holomorphic matrices Bj(£), v = —1,3, 
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j = 5,4,3, in an appropriate way. Omitting the straightforward but 
tedious computations, we present the result: 

V- |, V- hV- ^ ^n + mr, n,meZ: 




1 iV (oo) , 

r (°°) 1 _|_ 1 g(QO) r (00) j ' 

e 2"VWe(v + i±i)e(y+i)e(|)e(o) 
o(y+i)e(V)e(§)8'(i±^) 



y = i : iL. a (e) = r^( 



(oo) _ (do) = 2e(v)e(v+^)e'(i±^)e(o) 
q " w„e(v+i)e(v+5)e(i)9(§) ' 



y = i : i2_ 1 (^) = r« a8 . (6.40) 

The relevant correction function satisfies the estimates 

v (0 -/4 O(r7/24) - l5l<C ° n8t ' (6 41) 

in the domain of the parameter xq = xt~ 3 / 2 described by 

\V - ^ -n-mr\ > Ct~n +e , n, m G Z, e = const > 0. (6.42) 

6.6. Large t asymptotics of the RH problem [T](z/ = 3). In con- 
trast, the computation of the left rational multiplier for v = 3 
yields 

R 3 (0 = tl^I, (6.43) 
together with the additional condition, 

e(V) = 0, (6.44) 

or, equivalently, 

V = + n + mr . (6.45) 
In this case, the "external" parametrix is elementary, 

sf^Cfl = (ft(0r^(*s + vi)J«®", (6.46) 

and the correction function x(0 satisfies the estimates (I6.4ip . 

Remark 6.2. The above computation shows that, as t — » +oo, for each 
value of the deformation parameter x, xt~ 3 ^ 2 e _D ; the RH problem [T] 
is solvable either with v = — 1 or v = 3. Thus there is no room for 
the solvability of the RH problem with v = —5. Therefore, at least in 
the large t limit, the special solution to Pj has no triple pole collisions 

(2 1) 

and the corresponding solutions to PJ have no branch points. The 
relevant Malgrange divisor consists of smooth branches only. 
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7. Pole asymptotic distribution 

To compute the large t asymptotic distribution of poles of the special 
solution to Pj, we use the phase shift given in ( 16. 45ft and the definition 
of Fin flP7D , 

^-t 7/ \rA - B) = n + i + (m + \)r. (7.1) 

Using the canonical dissection of the Riemann surface, the difference 
UbA — u a B is expressed in terms of a single contour integral over the 
contour C depicted in Figure I7.1[ 



Figure 7.1. The contour for the canonical dissection of 
the elliptic curve Y. 

cubA — co a B = u a <j) U(z)fj,(z) dz. 

Inflating the contour C, we transform it to a contour encircling the 
infinite point which is a branch point of the curve. Then expanding 
the integrand at the infinity and using the residue theorem, we find 

tA-B = j>{U{z) - U^niz) dz = -£(3:ro + f£i), (7.2) 

and the pole distribution formula (17. ip yields 

3x + f 6 = K 7/4 {{n + |H + (m + |)w 6 ) . (7.3) 
Since £i is determined by xo via the Boutroux equations ( 16. 8p . equation 
(17.31) determines the pole position, x^' m ^ = x n)m t~ 3 ^ 2 as a transcendent 

function of two integers (n, m). As t — > +oo, the particular pole x "' m ^ 
approaches the attracting point satisfying the equation 

Zi(x ) = -|x . (7.4) 

— (2 1) 

7.1. Quasi-stationary solutions of equation PJ ' . Although it is 
not known how to solve the transcendent equation ( 17.41) . the problem 
of finding the attractor to the pole distribution significantly simplifies 
observing that the pole attractors correspond to the quasi-stationary 
solutions of equation Pj 2 ' 1 ' 1 , 

a(t) = o 0O * 3/2 (l + O(r e )), e>0. (7.5) 

Substituting (17. 5p into Pj ' ( 12. 7p . we find an algebraic equation for 
the parameter a^, 

, 

243 oo 1 2187 



4 _ 236 2 , m. = C7 61 

oo 243 oo ~ 2187 ' v'- u / 
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with the roots G {±|^=, i 2 -^}. However only one of the roots, 

a. = gf , (7.7) 

is consistent with the described above properties of the large t asymp- 
totic spectral curve for the special solution of P 2 . The linearization 
of equation PJ ' at the O-parameter power series solution with the 
leading order term coefficient (17. 7p has four linearly independent so- 
lutions. Two of them are exponential, ~ exp[±z^^(|) 3//4 |t 7//4 ], and we 
set them aside. Two other solutions of the linearized equation, ~ t -1 / 4 
and ~ t -1 / 4 In£ are relevant to our quasi- stationary behavior of the 
poles. Using them, we form the 2-parameter series 



a(t) = t 3/2 ]T(H \n(r 7 *)) k a k (t), a k (t) = 



k=0 



2VE 

9^3' 

fliOiOoi £ C are arbitrary, 

a 20 = — ^f^lOi a n = ~~ ^39?( a oi ~~ 2oio)aioj 

3v / 3(v / 15-50ag 1 +200a iai +3130af ) 
002 ~ 781575 ' 

° 30 = 307328 a 10' a21 = 307328 (1630qi + 861aio)a 2 Q, 

au = 98^9e(1907^ + 5216< + 55104a 01 ai + 233392a 2 )a 10 , 
a 03 = 3073180 ( 1630a oi + 25830e&ai + 218805a ia 2 + 788362a? ) 

+ ^rlSvs( 28605a oi + 258143a 1 o), 

(7.8) 

Let us relate the free parameters a i , aio in (17. 8p with the integers n 
and m in (17. 3p . Recall that, along the boundary dD + , two branch 
points £4 and £5 of the asymptotic spectral curve coalesce. Namely, if 
t — > +00, the limiting values corresponding to the attracting point are 

The asymptotics of the branch points in the model elliptic spectral 
curve compatible with the expansion (17.81) is given by 



£x = £2 = + |(ooi + 7a 10 + a 10 ln(r 7 / 4 ))f 7 / 4 + 0(r 7 / 2 In 2 t), 
£3 = - f (aoi + 4a 10 + a 10 ln(r 7 / 4 ))r 7 / 4 + 0(r 7 / 2 \n 2 1), 

- !(a 01 + How + ai ln(r 7 / 4 ))r 7 / 4 + G(r 21 / 8 ln 2 t). (7.10) 
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The asymptotics of the periods u a ^ as £ 5 — £ 4 — > have the forms 

". = ^(i-^+0((fc-ei) 2 )), 

u b = ^ In TB |^ I + 0((& - &) ln(& - 60). 6 - & -> 0. 

(7-11) 

Thus, for large £, 



w -^( 1 + r 7/8 )). 

wfc = -^lnt + ^ln^^ + 0(r^liit) > t^+oo. 

(7.12) 



Finally, the coefficients a±o and aoi determining the asymptotic series 
for a(t) (J7S1), 

a(t) = fjt 3 / 2 + rV* (a 01 + a 10 ln(r 7 / 4 )) + . . . , 

follow using the asymptotic formula (17. 3p . 



a 10 ln(r 7 / 4 )+a 01 + ---= (-|^-|6)t 7/4 +^((n+iK+(m+|K) = 
— — |oioln(t 7 ^ 4 ) — |a i — |aio 

+ 4(n + + i(m + i)^ N*-"*) + 1° 3$*] + • • ■ 

(7.13) 

Equating coefficients of ln(t~ 7//4 ), we find aio = ai™' 7 ^, 

a ft."> = ( m+ I)_^_ mGZ + , (7.14) 

while the constant terms yield the family of values of aoi = a^ ,n \ 

m G Z + , n G Z. (7.15) 

The formulas (17.141) and (I7.15P with (17. 8p yield the asymptotic formula 
(11.81) for a^ 171 '^ (t) which implies that the poles of the special solution to 
Pj in a vicinity of the attracting point x* = |^| form a regular lattice 
with the slowly modulated intervals and a boundary formed by the line 
of poles corresponding to the values m = and n G Z. In particular, 
the interval between two the most right vertical lines of poles is given 
by 

a (o,n )(t) _ a (i,„, (t) = _ r i/4^ ln[r 7 / 4 ^^] + Oit-^ In 2 1), 

Note that the boundary dD + formally corresponds to m = — |, n G 
R, and the distance between the first vertical line of poles and dD + is 

a*(f) - a<°">(f) = -t-^^H^^T^] + 0{rVt). 
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8. Problems and perspectives 

(2 1) 

Above, we have presented and studied equation Pf , to our best 
knowledge, the first of the differential equations that controls the isomon- 
odromy deformations of a linear ODE with rational coefficients and 
does not possess the Painleve property. For the first glance, its exis- 
tence contradicts the theorem by Miwa and Malgrange. However, it is 
not the case. The absence of the Painleve property in Pj 2,1 ' 1 is related 
to the fact that the domain of the solvability of the corresponding RH 
problem in the 2-dimensional complex space with the coordinates (t, x) 
is restricted to the Malgrange divisor of Pj , i.e. to the set of the complex 
lines (t, a(t)) which are allowed to intersect. Actually, this fact provides 
us with the important information on the nontrivial analytic structure 
of the Malgrange divisor for Pj which forms a Riemann surface with 
the infinite number of sheets and third order branch points. 

(2 1) 

Besides this, the discovery of equation Pj ' provides us with a 
new wide field of research. For instance, it is interesting to explore 
the possibility of existence of similar equations associated with other 
isomonodromic solutions of KdV or other integrable PDEs like Nonlin- 
ear Schrodinger or Pohlmeyer-Lund-Regge equations. 

Other interesting problem not discussed above is the structure in 

(2 1) 

Pj ' induced by the singularity reduction of the Hamiltonian structure 
in Pj 2 (and, in the case of a successful extension of the singularity re- 
duction methodology to the hierarchies associated with other Painleve 
equations, the structures induced by the Weyl symmetries). 

Finally, we mention the problem of characterization of such initial 
data to Pj (and other isomonodromy deformation equations) whose 
singularity reductions do not have branch points. In the Pj 2 case, we 
conjecture that its special solution considered above does not have 

(2 1) 

merging poles for any t and therefore the relevant solution of Pj 
does not have branch points. If this conjecture is true, it can serve as 
one more characteristic property of this special solution. 
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